Elementary Excitations in Quantum Antiferromagnetic Chains: Dyons, Spinous and 

Breathers. 
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Considering experimental results obtained on three prototype compounds, TMMC, CsCoCla (or 
CsCoBrs) and Cu Benzoate, we discuss the importance of non- linear excitations in the physics of 
quantum (and classical) antiferromagnetic spin chains. 
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I. INTRODUCTION: 

In the last twenty years, various experimental and the- 
oretical studies have been devoted to antiferromagnetic 
(AF) chains. One fascinating problem concerns the na- 
ture of the elementary excitations in such systems. As 
a crucial ingredient, one has to refer to the concept of 
non-linear excitations (NLE). Typical examples of NLE 
are provided by the solitons and the breathers. These 
particular excitations are solutions of the sine- Gordon 
equation, which, in the wide field of the NLE, is to 
be considered as a prototype model The first ev- 
idence of fluctuations associated with solitons in AF 
chains was obtained on a classical spin system (in the 
80's) iQ. Recently, evidence for breathers was experi- 
mentally obtained on a quantum spin system |^ . A short 
review is presented, along which, the role of the NLE in 
AF chains is discussed. Classical as well as quantum 
spin systems are analyzed. The case of anisotropic spin 
models is considered and a relation with the isotropic 
case is proposed. Quantization effects are also shown 
to complete the descriptions of the NLE in AF chains: 
for the solitons, this leads to another NLE concept, 
the dyons; for the breathers to a discretization of the 
excitation spectrum. Along this paper, we shall refer 
explicitly to three compounds: (CH3)4NMnCl3, alias 
TMMC, two Co compounds, CsCoClg and CsCoBrg [|, 
and Cu(C6H5COO)2.3H20, alias Cu Benzoate §. All 
of them are known for very long to be good models of 
1-dimensional antiferromagnetic chains. In TMMC, the 
spin value is large, s = 5/2, and a semi-classical spin de- 
scription is a good approach. The other compounds, how- 
ever, are representative of quantum spin systems, with 
s = 1/2. 



II. TMMC: 

TMMC is known to provide examples of quasi- isotropic 
Heisenberg AF chains. At Tco — 20 K, however, because 
of the small dipolar anisotropy D, a crossover occurs, 
which changes the spin system from isotropic to easy- 
plane behaviors The spins can be viewed as being 
repelled in the XY plane perpendicular to the chain Z 
direction. Then, if one applies a magnetic field H in 
the XY plane - along the Y axis, for instance - another 
crossover occurs because of the preferred spin- flop config- 
uration. In this later phase, the spins are mainly aligned 
along the X axis: they form an Ising-like AF state where 
the soliton regime takes place. The soliton phase diagram 
for TMMC is shown in Fig. 1 [§. 
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FIG. 1. Soliton phase diagram of TMMC [4]. Hc^ is the 
cross-over field between the soliton regimes A and B (see the 
insets) . 



1 



Two soliton regimes (A and B in Fig. 1, see also the 
insets) can be defined. Both they correspond to the sine- 
Gordon model. The difference between A and B results 
from the relative "competition" between the two effective 
anisotropics, one due to D, the other induced hy H 
We mainly limit the discussion to the "transverse" soliton 
regime A. For TMMC, the initial Hamiltonian is given 
by: 

H = Y, 2 Js„s„+i - 2I?«+i - gfiBHsl (1) 

n 

with J ~ 6.8 K and D ~ 0.16 K and s = 5/2. Each spin 
can be described as a vector, the orientation of which is 
defined by two angles, 9 and </). In configuration A, the 
spins are assumed to remain within the XY plane. Ac- 
cordingly, remains fixed to the value 7r/2, and 4> is de- 
fined by the angle between the spin and the external field 
H. In the Ising-like phase of TMMC, the spin vectors is 
simply written as — (— l)"cos(0). Within a classical 
description (and in the continuum limit), and after the 
variable change = — 2(/) is made, hamiltonian (1) is 
changed into 

/+00 
dz/2 
-OC 

[l/C^{dil!/dtf + {d^/dzf + m2(l - cos(V'))] (2) 

which is that of the sine- Gordon model. Here, C = 
4Js[l — Z?/(2J)]^/^ is the soliton velocity and m = 
gfj^sH/AJs defines the soliton mass. The elementary ex- 
citations of this model are typical non-linear excitations 
(NLE): the solitons (and antisolitons) and the breathers. 
A soliton in TMMC (or an antisoliton) can be viewed as 
a domain wall (extending over several lattice spacings), 
with a well-defined shape as represented in Fig. 2. The 
breather modes are soliton- antisoliton bound-states. As 
a general property, when NLE move along the chains, 
they maintain an exact balance between two energy con- 
tributions, the "potential" and the "internal" energies 
defined by the two last terms of (2), respectively. This 
explains the remarkable "integrity" of the NLE excita- 
tions ^. 




FIG. 2. Schematic representation of a (broad) soliton in 
TMMC. 

As can be seen in Fig. 2, inside a soliton, the spins 
undergo a tt rotation (this provides an alternative and 
a more general definition of solitons in AF chains). By 



looking at the fluctuations induced by the soliton mo- 
tion [HJ^, an accurate description of the soliton regimes 
in TMMC has been obtained. Concerning the breather 
excitations, the situation is experimentally more compli- 
cated. In spin systems, a discretization of the breather 
spectrum is to be expected, which, then, should result in 
a set of discrete excitation branches, hereafter denoted Bi 
... B„. In this quantization process, the lowest breather 
Bi coincides with the usual magnon excitation Q. If we 
consider only the A soliton regime, the lowest breather 
gap is that of the lowest magnon branch. The experimen- 
tal gap value is in agreement with the value predicted 
by (2): = gusH. In TMMC, however, there exist 
two competitive anisotropics as mentioned above. Ac- 
cordingly, there is another magnon branch with another 
gap: = AS^/D/J 0.8 meV). This higher gap is 
that of the lowest breather branch associated with the B 
soliton regime (see the insets in Fig. 1). In the quan- 
tization process, the breathers must be considered as a 
whole. As a result, additional peaks can be expected in 
the breather spectrum. This agrees with the observation 
of "unexpected" peaks at energies E ~ E^ ± Ef" ^ . 



Villain model : 

a) Groundstate : + - + - + - + - + - + - + - + - + - + - + - + 



c) Dyon : 




FIG. 3. Groundstate and elementary excitations in the Vil- 
lain model (hamiltonian (3)). 



III. THE CO COMPOUNDS: CSCOCL3 AND 
CSCOBR3: 

An alternative soliton (or kink) model for AF chains 
has been proposed by J. Villain, in 1975 ^. It apphes to 
the following quantum (s = 1/2) Ising-like Hamiltonian 

H = 2J5] + £(«+! + slsl^,) (3) 

n 

with e << 1. For this simple model, the groundstate is 
the Neel state as represented in Fig. 3a. 

The first excited state is displayed in Fig. 3b: it con- 
tains a narrow domain wall (extending over only one lat- 
tice spacing) and, as expected, it is associated with a tt 



2 



rotation of the sublattices. For this model, the dispersion 
of the soliton branch is shown in Fig. 4a. 
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FIG. 4. a) Dispersions of the soliton excitations in the Vil- 
lain model: the solid line is for a small e value; the dashed 
line, for a large e 1). The arrows illustrate experimen- 
tal transitions (see text), b) The fluctuation spectrum of the 
(individual) soliton mode (in zero field). 

In the same figure, the arrows show the transitions to 
be realized in experiments. The transitions induced from 
the groundstate (the dash arrow) are forbidden. They 
would require that an infinite number of spins be simul- 
taneously flipped (in order to realize the tt rotation of 
the sublattices). This establishes a pertinent result for 
Ising-like AF chains: the dispersion of a soliton branch 
cannot he determined experimentally. Only transitions 
inside the soliton branch (shown by the full arrow) are 
possible: they describe the fluctuations of (individual) 
solitons when they move along the chains. The corre- 
sponding spectrum is drawn in Fig. 4b: the fluctuations 
are seen to develop at very low energy (around zero en- 
ergy). A peak (i.e., a single "maximum") in the fluctua- 
tion spectrum is expected at Emax{q) = ^eJs\n{q) |^,|ll|. 




FIG. 5. Zeeman splitting of the soliton dispersion in the 
Villain model for a field parallel (a) and transverse (b) to the 
Ising direction. 

Such a kink in quantum spin chains is, in fact, a degen- 
erate state. As shown in Fig. 3c, it can be defined in two 
ways, which both agree with the tt rotation defined above. 
They differ, however, by the spin value (±1/2) inside the 
domain wall. A kink, i.e., a soliton, is a doublet state, 



entirely defined by the additional quantum spin number 
5 = 1/2. Accordingly, in a field a Zeeman splitting 
occurs in the dispersion, as shown in Figs. 5a and b. 
This splitting is different for H applied parallel (iJ|| ) and 
perpendicular (iJj.) to the Ising axis Z. The allowed tran- 
sitions are now induced between the two split branches 
(the full arrows in Figs. 5). Finally, this splitting results 
in a "doubling^^ of the observable soliton modes with two 



maxima, as shown in Figs. 6a and b |I2 
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FIG. 6. Doubling of the soliton modes in the Villain model, 
for H parallel (a) and perpendicular (b) to the Ising direction. 

The Co compounds, CsCoCla and CsCoBra, are good 
examples of the Villain model, with J ~ 75 K and 
£ ~ 0.12. The fluctuations associated with the solitons 
have been clearly identifled. In particular, with CsCoBrs, 
the low-energy fluctuations characterizing the (individ- 
ual) soliton modes have been observed by neutron in- 
elastic scattering (NIS) jl3|. Evidence for a "maximum" 
in the fluctuation spectrum requires a high instrumental 
resolution. Measurements have been performed at the In- 
stitut Laue-Langevin (ILL, Grenoble, France) using the 
high-flux high-resolution spectrometer IN14 ||lj]. An ex- 
ample of the soliton mode observed in zero field for the 
wavevector transfer Afc = q = tt/2 is displayed in Fig. 
7a. A single "maximum" is clearly detected. Its position 
agrees well with the expectation (the full line is a theo- 
retical curve which takes into account the instrumental 
resolution, 6E ~ I meV). In Fig. 7b, the same energy 
scan performed in presence of a magnetic field {H ~ 9.82 
T) is reported. Two maxima are now visible: this re- 
sult (the full line is the theory) establishes the predicted 
doubling of the soliton modes. 

The neutron data reported in Fig. 7 correspond to 
magnetic transitions associated with the particular mo- 
mentum transfer Afc ^ q = 7r/2 (the full arrows in Figs. 
4a and 5). Transitions with Afc = q = (the vertical 
dash arrows in Fig. 5) can also be considered. They are 
realized in an ESR experiment. Such measurements have 
been performed on the compound CsCoCls ||l^. They 
are referred to as Soliton Magnetic Resonance (SMR). 
Examples of SMR signals are reported in Fig. 8, for 
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both and H±. As shown by the full lines (theoreti- 
cal curves), the SMR data agree also very well with this 
concept ol & S — 1/2 soliton state. 
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FIG. 7. NIS on CsCoBrs: a) Soliton mode in zero field; 
b) doubling of the soliton mode in a transverse field. These 
energy scans correspond to the wavevector transfer q = -k 12. 

From these NIS and SMR measurements, an important 
property is established: in quantum Ising-like chains, a 
soliton is a S* = 1/2 doublet excitation. In fact, this result 
generalizes to any larger spin values, i.e., s > 1/2. This 
has been established, first, by Haldane in a semi-classical 
approach |l6): a soliton in an AF chain contains inter- 
nal degrees of freedom. After quantization, they result 
in a discretization of the soliton states. At this point, a 
drastic difference is established between integer and half- 
integer spin chains: in the former case, the lowest soliton 
state is a S" = excitation (this is the basis of the "Hal- 
dane conjecture"), while in the later case, it is a = 1/2 
excitation (in full agreement with the results on the Co 
Compounds). 




FIG. 8. SMR signals observed in CsCoCla, for H parallel 
and H perpendicular to the Ising direction. 

Hence, in half-integer spin systems, because of these 
additional degrees of freedom, the domain walls repre- 
sented in Figs. 2 and 3b, are not only solitons: they 
are dyons, according to Affleck ||l^. The dyon concept 
is used in particle physics: it describes a particle with 
an additional "charge" (cf: the magnetic monopoles). In 
the present NLE context, the S" = 1/2 solitons observed 
in the Co compounds are dyon analogs in quantum AF 
chains 0. 



IV. QUANTUM ISOTROPIC SPIN CHAINS : 

Up to now, we have considered Ising-likc chains, i.e., 
anisotropic systems, where the NLE are rigorously de- 
fined. Within the same context, we may try to reach a 
physical picture for isotropic chains. We follow the dis- 
cussion presented by Haldane for classical spin systems 
|l6t , but we apply it to the quantum case of the Villain 
model. Let us consider the change occurring in the soli- 
ton dispersion while the anisotropy is slowly removed, 
i.e., when e — > 1 in Eq. 3. As shown by the dash line in 
Fig. 4a, the amplitude of the dispersion increases, and 
at £ = 1, i.e., the isotropic point, the gap of the soliton 
branch would close (at the momentum value fc = 7''/2). 
A phase transition occurs ||T^, which can also be viewed 
as an isolant-metal transition ||l^. The dyon states (i.e., 
the pairs of -1-1/2 and —1/2 solitons) are seen to fill in 
the groundstate of the new phase, which is that of the 
quantum isotropic Heisenberg hamiltonian 

B^Y^ 2 Js„s„+i (4) 

n 

For this model, another concept, the spinon concept, is 
commonly used [|2^. A spinon is defined as an "entity" 
associated with a ±1/2 spin value and the groundstate 
of this isotropic hamiltonian is made of pairs of spinous. 
Such pairs of spinous compare well with the ±1/2 dyon 
states defined above. 

The description of the properties of quantum Heisen- 
berg chains started with the Bethe ansatz, in 1931 pl| . 
In the 60 's, based on the Jordan- Wigner transforma- 
tion, many studies were referring to a model of strongly- 
interacting fermions In the 70's, more sophisticated 
analyzes were proposed [ p3[ , which yields the modern 
field-theoretic derivations, which in the recent years, have 
renewed and enlarged our understandings of so many 
quantum spin systems. Cu Benzoate is such a recent 
example. For this compound, one is led to refer also to 
the sine-Gordon model [p4| . 

In anisotropic systems, the spin direction is easily de- 
fined (for instance, by the angle ^ in Fig. 2), and/or 
with respect to the groundstate. The pertinent variable 
in the sine-Gordon hamiltonian of Eq. 2 remains this 
angle (/>. In isotropic systems, however, there is no defi- 
nite direction for the spins. Alternative descriptions have 
been proposed |2^ . Usually, they starts with the Jordan- 
Wigner transformation. This leads to a fermionic repre- 
sentation of the hamiltonian with, in reciprocal space, the 
fermions operators a, (a+) and /3g (/3^). It is then recog- 
nized that the significant physical quantities are not con- 
tained in the fermion operators but rather in the fermion 
densities: p\ ~ J^k <^k+q^k and p' ~ J2k Pt+q(^k (this 
is the so-called "bozonization" ) . Another transformation 
- this is the basis of the field theoretic approach - brings 
back the equations to the real space where a field "^/i" is 
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defined. Applying this procedure to the isotropic hamil- 
tonian (4) leads to 



H ^ J/8 



+ 00 



dzl/2 l/C'^{di>/dtf + {dip/dzf (5) 



with C = ttJ. In this expression, the field ij} does not 
compare to the angle if) introduced in Eq. 2: ip is re- 
lated to the fermion densities p+ and defined above. 
Eq. 5 is not the sine-Gordon Hamiltonian. The last 
term which, in Eq. 2, introduces the non-linearity is 
missing. Accordingly, the excitations corresponding to 
Hamiltonian (5) follow /mear dispersions (i.e., no energy 
gap). They describe the lowest part of the spinon contin- 
uum, which characterizes the excitation spectrum of (4). 
As it is known, the application of a field H on hamilto- 
nian (4) develops a dynamical incommensurability'^: the 
wavevectors where the critical zero-energy fluctuations 
take place are shifted by Sqinc = 2iTa, where a, which is 
the magnetization per spin, is an increasing function of 
H. The dispersions, however, remain linear at low energy 
and no energy gap is induced by a field. 
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V. CU BENZOATE: 

In the 70's, CuBenz was considered as a good example 
of quantum isotropic AF chains p5| ]. Recently, however, 
it has been observed that, in a field, an energy gap opens 
in the energy spectrum in contradiction with the 
expectation for the isotropic model. According to Affleck 
and Oshikawa |^^, this behavior is due to the presence 
of an additional "staggered" field, hst [0. In a field, the 
effective hamiltonian for Cu Benzoate is to be written: 



H 



E 



2 Js„s„+i - gp.BHsl - gp^shstSn 



(6) 



with J ~ 5.8 K, and where h^t is proportional to the 
applied field: hst — 0.09iJ (for H applied parallel to 
the c axis) |2^] . The field-theoretic derivation of this 
hamiltonian [ P4[ transforms Eq. 6 into the sine-Gordon 
equation: 

/+00 
dz/2 
-00 

l/C^id^j/dtf + [di^/dzf + m^{l - cos(V^))j (7) 

where C and the field are defined in a similar way as 
in Eq. 5. The soliton mass is now a function of the 
staggered field: fh :^ l.85{hst/ J)^/^ \Hhst/ J)\ H. The 
low-energy excitations of (7) are solitons (antisolitons) 
and breathers. 



FIG. 9. A model for the dispersion of the elementary exci- 
tations in Cu Benzoate: the soliton/antisoliton branch gives 
a gap a.t q — and at the incommensurate wavevector 
q — n — Sqinc ■ The breather branches give a gap at q = n and 
at the incommensurate wavevector q = Sqinc- The dashed 
arrows corresponds to NIS [26], the full arrows to ESR [3]. 

In Fig. 9, a complete representation of the elementary 
excitations in Cu Benzoate is proposed. The bold lines 
correspond to the sine-Gordon predictions (low-energy 
limit). In this figure, three breather branches are repre- 
sented: Bi, B2 and B3. 

We now review briefly the recent experiments per- 
formed in Cu Benzoate. In Fig. 9, the transitions ob- 
served by NIS are shown by the dash arrows (these mea- 
surements have been performed at a single field value 
_ff ~ 7 T) [ p6[ . The dynamical incommensurability has 
been observed at q = n— Sqinc- At this precise wavevector 
value, the gap of the soliton/antisoliton branch is directly 
measured. At q — tt, the two first breather modes, Bi 
and B2, have been observed |^^. Recently, ESR mea- 
surements have also been performed in Cu Benzoate 
In fig. 9, the transitions observed by ESR are shown 
by the vertical dotted and solid arrows. The solid ar- 
row shows that the soliton/antisoliton gap is measured 
at 5 = 0. The dotted arrow probes the breather modes, 
also at g = . All these ESR measurements have been 
performed in a wide field range, up to 7? ~ 20 T. The ob- 
served field dependence of the soliton gap is displayed in 
Fig. 10. In low field, the agreement with the sine-Gordon 
model (the full line) is very good: Eq — fhJ ~ H^/^. An 
appreciable deviation, however, occurs in high fields. In 
Fig. 11, the field dependence of the three first breathers, 
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Bi, B2 and B3 as probed by ESR 
tained in the full field range. 



(i.e., at q = 0) is ob- 
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FIG. 10. Field dependence of the soliton gap observed by 
ESR (Voigt configuration) in Cu Benzoate. 



VI. CONCLUSION: 

With these compounds - i.e., TMMC, the Co com- 
pounds and Cu Benzoate - the importance of the NLE 
in the low-T properties of AF chains is well exemplified. 
They are seen to play a crucial role in both anisotropic 
and isotropic systems. They come in the description of 
both quantum and classical spin chains. TMMC and 
Cu Benzoate are both well described by the sine-Gordon 
model. Fundamental differences, however, occur between 
the two systems, in particular for the experiments. 




In Ising-like chains (TMMC and the Co compounds), 
transitions from the groundstate to the soliton state are 
forbidden (see the dash line in Fig. 4a). In such sys- 
tems, the presence of solitons is detected by looking at 
the low-energy fluctuations induced by the soliton mo- 
tion P,p^. In Cu Benzoate, however, the soliton disper- 
sion can be determined experimentally: the transitions 
from the groundstate to both the soliton and the breather 
branches are allowed, (see the arrows in Fig. 9). With 
Cu Benzoate, the discretization of the breather spectrum 
in a sine-Gordon model is now well established. 
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